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IA.BSTRACT: Most models for soft interactions which were proposed prior to the measurements at the LHC, 
^-H are only marginally compatible with LHC data, our GLM model has the same deficiency. In this paper we 
^ investigate possible causes of the problem, by considering separate fits to the high energy [W > 500 GeV), 
_^ and low energy {W < 500 GeV) data. Our new results are moderately higher than our previous predictions. 
^) Our results for total and elastic cross sections are systematically lower that the recent Totem and Alice 

published values, while our results for the inelastic and forward slope agree with the data. If with additional 



L^ experimental data, the errors are reduced, while the central cross section values remain unchanged, we will 
l^ peed to reconsider the physics on which our model is built. 
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1. Introduction 

The new LHC data on soft interaction scattering at higii energy (see Refs. [1-4]) is only marginally com- 
patible with updated Pomeron models [5-8], which have been successful in reproducing the cross section 
data in the lower energy range. The implication is that, as it stands, our understanding of long distance 
physics at very high energies is limited. From an optimistic point of view this may imply no more than 
the need to adjust the Pomeron model's parameters. From a pessimistic point of view, this may suggest 
the need for a comprehensive revision of the main ingredients of Pomeron models applied to high energy 
soft interactions. 

Specifically, the model we propose [5] is built using an input Pomeron with a relatively large fitted 
intercept Ajp = ajp — 1 = 0.2 and exceedingly small slope a'p = 0.02GeV~^. These values are in 
accord with AdS-CFT correspondence [9-13]. Note that in N=4 SYM [9], Ajp = 1 - 2/^/A ^ 0.1 - 0.3, 
corresponding to experimental estimates based on multiparticle production, as well as HERA DIS data [14] 
in which A ~ 5 — 9. The other basic ingredients of our model are the large Good-Walker(GW) [15] 
contribution to diffraction production, and a small Pomeron self interaction. Both are direct consequences 
of the AdS-CFT correspondence. 

If the present central values of the LHC cross section data points, which have relatively large errors, do 
not change significantly with the forthcoming, better statistics, measurements. This would suggest either 
of two extreme options: 

1) The new LHC data do not support the main theoretical concepts of our model. This may stem from our 
reliance on Reggeon calculus and pQCD, which led to a single Pomeron model, rather than the traditional 



distinction between a soft and a pQCD hard Ponieron, and/or from our realization of AdS-CFT N=4 SYM 

ideas. 

2) Our procedure for adjusting parameters may be deficient, requiring a more sophisticated data analysis 

which may yield satisfactory results. We note that the fitted data base [5] contains no LHC data. Moreover, 

the low energy {W < bOOGeV) total, elastic and diffractive cross sections which constitute the major portion 

of the fitted data points have rather small errors. Consequently, our fitting procedure is not well balanced 

as the main contribution to our x^/d-o.f. stems from the low energy data. 

An alternate, and probably a more realistic option would be based on elements originating from the above 

two propositions. 

In this paper we check the second option. To this end we removed the low energy data and only 
fitted the high energy data (W > bOOGeV), including the available LHC soft cross section data points, so 
as to determine the Pomeron parameters. Having adjusted these parameters, we tuned the value of the 
Reggeon-proton vertex, which enabled us to obtain a smooth cross section behaviour through the ISR-LHC 
energy range. We hope that this exercise will clarify to what extent our model has intrinsic deficiencies, 
or do we just have a technical problem in the procedure for adjusting our free parameters. 

The results of this paper can be summarized as follows: First, we show that in spite of the fact that 
the values of the parameters, extracted from our current fitting, are different from our previous values, the 
overall picture remains unchanged. Second, our updated total and elastic cross sections are slightly lower 
than the published TOTEM values [4], but still within the relatively large experimental error bars. Should 
future LHC measurements confirm the present TOTEM values, we will need to revise our dynamic picture 
for soft scattering. 

2. Our Model 

The ingredients and formulae of our model have been published (see Ref. [5]). However, in order to produce 
a self contained presentation, we start with a brief overview of our formalism. 

As we have mentioned, one of our main input assumptions is the GW mechanism [15], which plays a 
significant role in the calculation of the eikonal shadowing corrections. To this end we took into account a 
two channel formalism in which we introduced two eigen wave functions, V'l and '(/'2) which diagonalize the 
2x2 interaction matrix T, 

Ai,k =< i'i V'fc|T|V'j' ipk' >= Ai,k Si,i' Sk,k'- (2.1) 

In this representation the two observed states are an hadron (a nucleon in our calculations), denoted 
by the wave function iph and a diffractive state ipD- Note that, we replaced the rich population of the 
diffractive Fock states by a single state with unknown mass. This representation provides a considerable 
simplification of our calculations at the price of not being able to calculate the mass dependence of GW 
diffraction production. The two observed states can be written in the form 

iph = a ipi + /3 ■ip2 , ipD = -131^1 + 0:^^2, (2.2) 



where, 0^ + 0^ = 1. Using Eq. (^.l]), we can rewrite the s-channel unitarity constraints in the form 

2ImAi,fc(s,6) = |Ai,fc(s,6)|2 + G^ffc(s,6), (2.3) 

where, G**], is the contribution of all non GW inelastic processes. 
In a general solution of Eq. (|2.3|) 



A;^(^,6) = ^(l-exp(-^^^)), (2.4) 

G^;^fc(s,6) = l-exp(-0,,fc(s,6)), (2.5) 

in which fij^fc are arbitrary. In the eikonal approximation Jlj^fc are real and amplitude Aik are pure imaginary. 
In general we have 4 Aij. amplitudes, however, for pp and pp Ai^2 = ^2,1 • From Eq. ( p.5[ ) we deduce that 
the probability that the initial state (i, k) remains intact during the interaction, is P^j^ = exp (— fij ^(s, h)). 

The input opacity Oj ^(s, b) corresponds to an exchange of a single bare Pomeron. 

^^,k{s,b) = gi{b)gk{b)P{s). (2.6) 

P{s) = s and gi{b) is the Pomeron-hadron vertex parameterized in the form: 

g^ (6) = gi Si{b) = ^mfbKi {rmb) . (2.7) 

47r 

Si{b) is the Fourier transform of , ^i — jw; where, q is the transverse momentum carried by the Pomeron. 
In our calculations we assume that the slope of the Pomeron trajectory is* a'jp = 0. This is compatible 
with the exceedingly small fitted value of a^, and in accordance with N=4 SYM [9]. 

In our model [5], the Pomeron's Green function that includes all enhanced diagrams is approximated 
using the MPSI procedure [16], in which a multi Pomeron interaction (taking into account only triple 
Pomeron vertices) is approximated by large Pomeron loops of rapidity size of Ins. 

We obtain 

Gr (K) = 1 - e.p (^) ^ r (0. ^) , (2.8) 

in which: 

T{Y) = 7 6^^"^. (2.9) 

r (0, l/T) is the incomplete gamma function (see formulae 8.35 in Ref. [17]). 

Summing the net diagrams [5], we replace gi{b) by a more complicated vertex function which, together 
with the enhanced diagrams, results in the following expression for i}i^k{s,b): 

,, f , m(b') gk(b-b') (1/7 Gjp(T(y))) 
n'^{Y;b) = cfb' ^^ ^ \\ ,. ./ .. . (2.10) 



1 + iG3iph)Gjp(T{Y) 



gi{b']+gk{b-b 



*Actually, we sum the diagrams for the Pomeron interaction considering a'p — 0, but introduce a'p for the Pomeron 
exchange. Since, the output of our fit gives a small value of a'p ~ 0.02 061^"^, we consider that this procedure is justified a 
posteriori. 



Gsjp is the triple Pomeron vertex, and 7^ = / -^;^G'l^]p- Note we consider 7 as an independent parameter 
in our fit to the data. 

For the elastic amplitude we have: 

a^iih) = (a^^i,i + 2a^p^Ai^2 + /3^^2,2) • (2.11) 

For diffraction production we introduce an additional contribution due to the Pomeron enhanced mechanism 
which is non GW. For single diffraction we have (see Fig. la): 

41/ = / d^b' 2 ^ (^ \) 9i {b - b', mi) gi (b', mi) g^ (b' , rrik] 

X Q(^gi,mi,b-b',Ymj Q (^gk,mk,b',Y -Y^j Q [gi,mub' ,Y -Y^ , (2.12) 

where, 

Q{g,mMY) = i^f^r\v^ e ^ . V (2-13) 

1 + {G3JP/V 9 Gf [Y) S [b, m) 

The structure of Eq. ( 2.12| ) can be understood from Fig. |T|-a. Q {g, m, b; Y) describes the sum of the 'fan' 



Pomeron diagrams. As shown in Fig. ffl-a, we have one cut Pomeron in Fig. ffl-a, which we express through 
the Pomeron without a cut, using the AGK cutting rules [18]. 

For double diffraction we have (see Fig. lb): 

4,1 = j dH'Agi{b-b\m^ gt (b', m^) 

X Q(g„mi,b-b',Y-Yi)e^^^''Q(gk,mk,b',Yi-5Y). (2.14) 



This equation is illustrated in Fig. |l|-b, which displays all ingredients of the equation. We express each 
of two cut Pomerons through the Pomeron without a cut, using the AGK cutting rules [18]. 

Eq. (2.12|) and Eq. (|2.14|) are the simplifications of the exact formulae of Ref. [5], which correspond to 



the diagrams of Fig. |l[ We checked that they approach the values of the exact formulae reasonably well, 
within 5 - 10%. 

For single diffraction, y = In (M'^/sqJ, where, M is the SD mass. For double diffraction, Y — Yi = 
In (M^/sq) and Yi — SY = ln(M2/so), where Mi and M2 are the masses of two bunches of hadrons 
produced in double diffraction, sq is the minimal produced mass, which is about 1 GeV . 

The integrated cross section of the SD channel is written as a sum of two terms: the GW term, which 
is equal to 



,CW _ I .2. 



/ (fb |a/3{-a2 Ai,i + (a^ _ /32) ^^ ^ + (3^ ^2,2} |^ (2.15) 

Ai^j^ are given by Eq. (|2.4| ). The second term describes diffraction production due to non GW mechanism: 



a,"j^W ^ 2JdY^jdH (2.16) 

{ a^ Al%^^ e-^^O-;") + a^p^Al%^2 ^-''^Ay^b) + 2 a^ /S^ Al^^^ ^-\(^^Ay^b)+n,,,{Y;b)) 




Q(Ym; •..) 





Q(Y-Y,;...) 



a; 




Figure 1: The set of the diagrams for single (Fig. |^a) and double (Fig. |i|-b) diffraction. The wave like lines denote 
Pomerons. The solid horizontal lines correspond to interacting hadrons. The vertical dashed line denotes the cut 
that corresponds to produced particles. 



The cross section of the double diffractive production is also a sum of the GW contribution, 

2 



GW 



2 o2 



d^h OL^ /3 



A 



1,1 



2^1.2 + A 



2,2 



to which we add the term which is determined by the non GW contribution, 



nGW 



d% 



\a'A\ 



dd^^-n,My;b) + 2a^l3^Af^e-''''^^^'''^ + /3^ Ai^e-""^'^^^' 



^"}^ 



(2.17) 



(2.18) 



In our model the GW sector can contribute to both low and high diffracted mass, as we do not know the 
value of the typical mass for this mechanism, on the other hand, the non GW sector contributes only to 
high mass diffraction. 

3. Results 



Using these formulae we fit the data with energies W = ^/s > 500 GeV, including the LHC data. Our 
best set of parameters is presented in Table 1. In the same table we also show the values of the parameters 
that we obtained in our previous pre LHC fit, which included the data for (W < 1800 GeV). Although, 
the two sets of parameters differ somewhat, the two key parameters, the Pomeron intercept Ajp and the 
slope of the Pomeron trajectory a'jp, remain almost the same. 

The comparison with the LHC data at W = 7 TeV is shown in Table 2. Even though we are definitely 
below the data points for total and elastic cross sections [4], the corresponding large experimental error bars 
do not enable a definitive assessment of our predictions. Our predictions for B^i and dinei are compatible 
with the TOTEM data. TOTEM's atot depends on the value of p = Re A/Im A, where Re A and Im A 
are the real and imaginary parts of the scattering amplitude A, and on daei/dt(t = 0). 

In our model we assume that Pomeron exchange leads to a pure imaginary amplitude. Since, we expect 
that the real part of the amplitude will be much smaller than the imaginary one, we can calculate the real 
part using a perturbative approach. 



Ajp 


/3 


a'jp{GeV-^) 


51 (Gey-i) 


52 (Gey-) 


mi ( GeV) 


m2 (GeV) 


0.21 (0.2) 


0.46 ( 0.388) 


0.028 (0.02) 


1.89 (2.53) 


61.99 (88.87) 


5.045 (2.648 ) 


1.71 ( 1.37) 


AjR 


7 


a'j^iGeV-^) 


9^{GeV-^) 


fff (GeF-i) 


Rli {GeV-^) 


G^p{GeV-^) 


-0.47 (-0.466) 


0.0045 ( 0.0033) 


0.4 (0.4) 


13.5 (14.5) 


800 (1343) 


4.0 (4.0) 


0.03 (0.0173) 



Table 1: Fitted parameters for our model. In parenthesis we include the values that we obtained in our model 
without LHC data (see Ref. [5]. The quality of the fit is x^ /d.o.f. = 0.86 (see the detailed explanation in the text.) 

First, we notice that for one Pomeron exchange 



Re^jp(s,6) 



tan ImAjp 



Having Eq. (3.1) in mind we can calculate the real part of the scattering amplitude as follows 



ReAjfc(s,6) = ReOjfe/2 exp 



Viik (s, h) 



(3.1) 



(3.2) 



We use Eq. (^T|) to calculate Rerij^. 

Our value of p is smaller than the COMPETE [23] value which was used by TOTEM. The COMPETE 
p fit is based on an extrapolation from data in the ISR-Tevatron range. The effect of our value of p being 
smaller than the COMPETE value, implies a change in atot of less than 1%. 

In Table 3 we present our results for atot, CTei, o'sd, o'dd, Bei, Bgd and amei- In parenthesis we put the 
values of these observables obtained from our previous fit. The new fit gives higher values for the single 
and double diffraction, while changes in all other observables are rather small. It is interesting to note 
that Block and Halzen [24] have "converted" [24] a recent measurement of the Pierre Auger Observatory 
collaboration of proton-air collisions with (Tf~'^^^ at W=57 it 6 TeV to ain for proton-proton collisions, and 
obtain the value of af^ = 90 ± 7{stat) ± 1.5{Glauber{+9/ - ll(sysi) mb and predict aff^ = 134.8 ± 1.5 mb 
at this energy. 

It is interesting to note that at W = 57 TeV, we have that the ratio -^^ = 0.74, while Block and 
Halzen [24] have the value 0.69, both far from the black disc value of 0.5. 

In Table 3 we show that we obtain the same value for the inelastic cross section, while for the total 
cross section we have 122 mb which is smaller than has been advocated in Ref. [24]. 

In Fig. m we show our present fit in the energy range 20 < W^ < ISOOGeV energy range. We see that 
the description is not as good as in our previous paper [5], nevertheless, our model describes the main 
features of low energy data as well. 



In Fig. y we plot the amplitudes Aik at different energies as a function of the impact parameter b. 
The structure of the b and s dependence, remains essentially the same as in our previous descriptions. 



w 


„raodel 
^tot 


exp 
^tot 


—model 
'^el 


<r 




7TeV 


94.2 mb 


TOTEM: 98.3 ±0.2*** ±2.8''2'"*mb 


22.9 mb 


TOTEM: 24.8±0.2"* ± 1.2"?'"*mb 




W 


—model 
^ in 


^Z' 


T:>model 
^el 


B:r 




7TeV 


71.7 mb 


CMS: 68.0±2"''''* ± 2.2'"™ ± 4^^*™p mb 
ATLAS: 69.4±2.4^^p ± 6.9^=^*™^ mb 

ALICE: 72.7 il.l™"'*^; _^ r^^extrap j^^3 

TOTEM: 73.5 ±0.6"* ± l.S''^''* mb 


19.8 GeV-^ 


TOTEM: 20.1±0.2"* ± 0.3"^"* GeV'^ 




W 


—model 
^sd 


exp 
'^sd 


—model 
'^dd 


exp 

'^dd 




7TeV 


IQ.b^^ + 2.6"^^ mb 


ALICE : 14.16 ± 3 mb 


5.98^^ + 1.166"*^^ mb 


ALICE: 8.86 ± 3 mb 





Table 2: Comparison of the predictions of our model with the experimental data at W= 7 TeV. 



y/s TeV 


1.8 


7 


14 


57 


atot mb 


75.6 (74.4) 


94.2 (91.3) 


104.0 (101) 


122.0 


(Tel mb 


18.2 (12.5) 


22.9 (23) 


26.1 (26.1) 


31.1 


asd{M < Mo) mb 




10.5 + (2.6)"'=^ (10.2) 


11.2 + (3.32)"G'^(10.8) 


12.8 + (3.91)"<^^ 


asd{M^ < 0.05s)mb 


8.97+ (1.95)"^^ (8.87) 


10.5 + (3.94)"'^^ (10.2) 


11.2 + (5.58)"^'^(10.8) 


12.8 + (8.19)"<=^ 


o'dd mb 


5.56 + (0.369)"^^ (4.46) 


5.98 + (1.166)"'=^ (6.46) 


6.55 + (1.5)"^^ (6.65) 


8.61 + (4.9)"^^ 


Bel GeV-^ 


17.6 (16.1) 


19.8 (19.3) 


21.2 (20.5) 


23.8 


B?r GeV-' 


6.36 


8.01 


8.78 


10.4 


(Tinel mb 


57.4 


71.7 


77.9 


90.9 



Table 3: Predictions of our model for different energies W. Mq is taken to be equal to 200Gey as ALICE measured 
the cross section of the diffraction production with this restriction. 



However, the amplitude An turns out to be smaller than that obtained previously [5]. This amplitude 
does not reach the unitarity bound even at W = 57 TeV. Recall that the unitarity black disc bound at a 
given (s, b) is only reached when Aii{s, b) = Ai2{s, h) = ^22(5, b) = 1. As seen in Fig. |^ this condition has 
not yet been satisfied at W = 57 TeV. In some sense, our model gives an example of why it is dangerous 
to apply the simple formulae that are valid for the black disc regime, for extracting information on the 
high energy scattering from the experimental data. In particular, as we have discussed at W = 57 TeV, 
we obtained atot = 122mb while a^i = 31mb in clear violation of the relation a^i = crtot/2 for a black disc. 

The impact parameter dependence shown in Fig. |3| reflects in the t dependence of the elastic cross 
section shown in Fig. |^. One can see that we are in agreement with the experimental results for both 
Tevatron and LHC energies in the forward cone for -t < 0.5 GeV^. In Fig. Q we show our prediction for 
W = UTeV. 

A significant check for our model will be the measurement of i-dependence for dagd/dt. As shown 
in Table 3 we predict quite a small slope Bg^ for the Good- Walker contribution of the single diffractive 
production. For large mass diffraction due to triple Pomeron interaction we expect that the slope will be 
about i3ez/2, since the slope for triple Pomeron vertex is rather small < I GeV^^ . For single diffraction we 
obtain a value of the slope Bgd = 8.01 GeV'"^ , (at t =0), this is compatible with the preliminary result of 



TOTEM as presented by Risto Orava 

As is well known, the survival probability for the large rapidity gap crucially depends on the b- 
dependence of the amplitude. To check our 6-dependence as well as for completeness of our presentation, 
we calculate the survival probability S'^ for a large rapidity gap using the general formulae of Ref . [26] . We 
obtain 

^2 = 9.76%(10%) at 1^ = l.STeV; S^ = 5.32%(6.3%) at W = 7TeV; S^ = 3.66%(4.4%) at 1^ = UTeV; 

(3.3) 
where we indicate in parentheses the values of the survival probability calculated in our previous paper [26] . 
Naively, we could expect that S"^ would be larger than in our previous model since the larger transparency 
of ^11 in this approach (see Fig. |3|). From Fig. ^ we see that Aii{b) in this model is smaller, and decreases 
faster than this amplitude in our previous model. Both these effects lead to an increase of the survival 
probability due to the Good- Walker mechanism. However, the considerable increase of the values for Gsjp 
and 7 lead to stronger screening, due to non Good- Walker mechanism, mostly due to contribution of the 
enhanced diagrams (see Refs. [5,26]). The second effect that leads to the decrease of 5^, is the value of a 
in this model, which is smaller than in our previous approach. Recall that a^ = 1-/32^ and 5^ oc a^, since 
the contribution of A12 and A22 to the value of the survival probability is negledgibly small due to large 
rii2 and r222 at small b. 

4. Conclusions 

We are able to reproduce the LHC experimental data within the present experimental errors and uncer- 
tainties. Therefore, we found a positive answer to the question formulated in the introduction: our model 
without any changes except the new set of parameters is able to describe the LHC data with a good ac- 
curacy. It is instructive to note, that we did not spoil the reasonable description of the low energy data 
(W < 500 Gey) being the only model on the market that is able to describe all available data on high 
energy scattering. We would like to stress that our predictions for atot and aei are below those published 
by TOTEM [4], while our values for amei and B^i are in agreement with the published experimental values. 
If the published errors will be reduced while the central values will remains approximately the same, our 
approach will need essential improvements and new ideas. 

As one can see from Table 1 the changes in parameters are not dramatic but three parameters: mi, 7 
and G3JP, turns out to be two to three times larger that in our previous approach [5]. These changes are 
driven by the LHC data and they deserve a discussion. The increase in 7 and G^jp is a direct consequence 
of the large diffractive cross section measured at the LHC. At first sight we could increase them even more 
to obtain a better description of the of the high mass diffraction production . Unfortunately, this is not 
true as the sum of enhanced and net diagrams depend strongly on the value of 7 and Gsjp, leading to a 
resulting decrease of the diffractive cross sections at larger values of these parameters. 

The most striking manifestation of larger value of mi we see in Fig. |3| which shows that at ultra high 
energy the amplitude An (b = 0) is still less that unity. The slow increase to the unitarity limit in An 



amplitude is a typical feature of our approach. However, the fact that this amplitude is less than 1 for 
6 = even at P^ = 57 GeV is certainly the consequences of the LHC data in the framework of our model. 

There are two lessons that we can learn from this fit. The first is that if the TOTEM collaboration 
will confirm their results for atot and a^i, we will need to reconsider the main physical ideas on which our 
approach is built. The second one, is the fact that one has to be very careful when considering the concept 
of the black disc regime. 
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Figure 2: Comparison with the experimental data the energy behaviour of the total (Fig. ||-a), elastic (Fig. ||-b), 
single diffraction (Fig. ||-c), double diffraction (Fig. ||-d) and inelastic (Fig. g-f) cross sections and elastic slope( 
Fig. ||-e) . The solid lines show our present fit . The data has been taken from Rcf. [f 9] for energies less than the 
LHC energy. At the LHC energy for total and elastic cross section we use TOTEM data [4] and for single and double 
diffraction cross sections are taken from Ref. [1]. 



11 



1 

0.8 




A, 


k 


W=1.8TeV 


r^ 




\ 


\\\2 \ 


0.6 


- 


^ 




\\ V\ ^22 


0.4 


"An 


\ 




V V 


0.2 


- 






>^^. . />^;^=^-.j... . /r>^ 



0.5 1 1.5 2 2.5 3 

bin fm 




Fig. §c 



Fig. §h 




1.2 



0.8 



0.6 



0.4 



0.2 



A, I 


W = 57 TeV 


:\ 


\a,2 \ 


■ \ 


\ \ A22 


■ A„ \ 


\ \ 


■ ^ 


V vV 



0.5 1 1.5 2 2.5 3 

b in fm 



Fig. |c 



Fig. Hd 



Figure 3: The impact parameters dependence of amplitudes Aik at different energies. The sohd lines are the results 
of this model, the dashed line are the amplitudes given by the model of Ref. [5]. 
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